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Abstract 

Nilsson et al. have recently shown how a velocity field in geostrophic 
and hydrostatic balance in an /-plane may be diagnosed from prescribed 
distributions of buoyancy and wind stress, in a basin with closed isobaths. 
I extend their analysis to cover basins with more complex depth contours, 
treating in particular the behavior of the flow in the presence of a saddle 
point. 

1 Introduction 

Recently, Nilsson et al. [2] showed that the steady circulation in a basin with closed 
depth contours, under conditions of geostrophic and hydrostatic balance, with arbi- 
trarily prescribed distributions of buoyancy and wind stress, may be determined by 
integral constraints, obtained by integrating around the contours. Taking account 
of bottom friction, they found that the net Ekman transport across the closed con- 
tours must vanish, and that this constraint determines the free barotropic motion 
and the entire velocity field. If, however, the depth contours are not all simple 
closed curves, but include one or more saddle points, the integral constraints may 
require the fluid velocity to be discontinuous on contours which include saddle 
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points. This discontinuity they showed to be inconsistent with their scaling as- 
sumption that the flow and the bathymetry vary over comparable length scales. 

In this paper we consider the steady flow in a basin which contains a saddle 
point. In the presence of linear bottom friction, the governing partial differential 
equation for the barotropic motion is shown to be elliptic and non-singular if the 
depth topography is smooth and has a bounded gradient. If we investigate the 
behavior of the solutions to the equation in the vicinity of a saddle point 

H(x,y) = H + Axy, 

we find that the equation is separable, and the solution may be expressed as a sum 
or integral of solutions to Hermite's equation. 



2 Mathematical development 

For a basin with variable depth H and with surface and bottom Ekman layers, 
containing water of variable density p r (l — q), Nilsson et al. [2 i derived the following 
equation for the steady, geostrophic current u in the interior: 

u(x, y, z) = 9 - R f V h q dz + ^-q b RVH + -^RVp , (1) 

/ J—H f fp r 

where g is the acceleration due to gravity, / is the Coriolis parameter, q b = 
q(x,y, —H), po is a barotropic pressure anomaly, and the linear operator R ro- 
tates vectors 90 degrees to the left in the xj/-plane. Nondimensionalizing so that 
p r = g = f = 1, we obtain, for the bottom velocity Ub (just outside the Ekman 
layer) : 

u b = q b RVH + RVpo- (2) 
Nilsson et al. [2] also show (using our non-dimensional notation) that 

u • VH = V • (u H) = -V • m, (3) 

where 

u = RVpo, (4) 
and the Ekman boundary-layer transport m is given by 

m = -R(r w - r b ), (5) 
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where r w is the wind stress and r b is the bottom stress, given, in a linear represen- 
tation, by 

n = h e u h , (6) 

where h e is the Ekman depth. 
From ((THBJ), we obtain 

V • (HKVpo) = V • [Rr w - h e R (q b RVH + RVp )} , 

which we may write as 

V • \(HR - h e I)Vp ] = V • (h e q b VH) -Vxr w , (7) 

I being the identity operator, or 

h e V 2 p + V# x Vp = V x r w - h e (q b V 2 H + Vq b ■ VH^j . (8) 

For non-zero h e , (JSJ) is elliptic. 

In the absence of friction (h e = 0), VH x Vpo = 0, and po may be an arbitrary, 
not necessarily continuous, function of H if the forcing terms on the right-hand 
side of dHJ are zero. The flow velocity is directed along the depth contours. 

For < h e H, if the depth contours are closed, Nilsson et al. [2] showed 
that uo is given by 

u = - [q(H) + f(H)] RVif, (9) 

where 

q=<[ q h \VH\ds(l \VH\ds) (10) 

JC(H) \JC(H) J 

and 



T(H) = h e - L i T w -ds[i \VH\ds\ , (11) 

JC(H) yc(H) J 

where integration is around the depth contour C(H). If the bottom contours 
are all smooth closed curves, with no saddle points, and t w is smooth, then p 
and the flow velocity will be smooth and continuous. If, however, the bottom 
profile has one or more saddle point, it is possible for (JHMll)) to give flow velocities 
which are discontinuous across a separatrix (the isobaths which touch at a saddle 
point). This violates the scaling assumption of Nilsson et al. j2] that the flow and 
bathymetry vary over comparable length scales, and it is likely that the solution 
to (jHJ) will have a boundary-layer structure in the vicinity of the separatrix. 
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Note that the zero bottom slope at a saddle point should not lead directly 
to any singularities in the solution, since if VH = 0, (jSJ) reduces to the Poisson 
equation 

h e V 2 p = V x r w . (12) 

3 Solution in the vicinity of a saddle point 

To determine the behavior of solutions to (jH)) in the presence of a saddle point, we 
consider a special case of the homogeneous version of with H = H^ + Axy, the 
saddle point being at the origin. We drop the zero suffix from p Q for convenience. 
We have 

h e V 2 p - Ax^- + Ay®?- = 0. (13) 
ox ay 

This linear partial differential equation is separable, with solutions of the form 
X(x)Y(y), satisfying versions of the Hermite equation 

^-X" - xX' + kX = 0; ^Y" + yY' - kY = 0. (14) 

There is thus the possibility of constructing boundary layer solutions by employing 
linear combinations of solutions to (|14j) for suitable combinations of values of k. 



4 Boundary- layer structure 



a. inner solution 

If we let x* = x/L and y* = y/L in (fT3HT4]) . where L = (2h e /A) 1/2 1 we obtain, 
removing the asterisks for convenience: 

V 2 p- 2*^ + 2^ = 0, (15) 
ox ay 

and 

d 2 X dX 

— T - 2x— + 2kX = 0; 16 

dx dx 

d 2 Y dY 

w + 2y __ 2kY = o. (17) 
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Equations ITBHT7I are variants of Hermite's equation, whose solutions for integer val- 
ues of k include Hermite polynomials and repeated integrals of the complementary 
error function pQ. 

For k > 0, Eq. El has solutions of the form 

X k = A x R k {x) + B{I k+1 e x \ 

Y k = A 2 l k erfc(y) + B 2 l k erfc(-y), (18) 

where A±, Bi, A 2 , and B 2 are constants, H k is the Hermite function of order k 
(reducing to the Hermite polynomial of degree k if k is an integer), 

erfc(z) = -= / e- z (19) 

\/7T Jz 



is the complementary error function, and X k denotes the following integral operator 
applied k times: 

POO 

lf(x)= / f(t)dt, (20) 



suitably generalized for non-integer k. For k < 0, the solutions are 

X k = C{I~ k eiic(-x) + £>iT - *erfc(a;), 
Y k = C 2 H_ k (-y) + D 2 l- k+1 e y 2 , (21) 

where C*i, D x , C 2 , and D 2 are constants. 

For the particular cases k = ±1 with B\ = D 2 = 0, and A\ = C 2 = 1, we 
obtain 

Xi = x, 

Yi = A 2 (-j/erfcfe) + (tt)- 1 / 2 ^ 2 ) + B 2 (yeric(-y) - (tt)- 1 / 2 ^ 2 ) ; 
X_i = d (xerfc(-:r) - (vr)^ 2 ^ 2 ) + A (-zerfc(x) + (Tr)~ 1/2 e - x2 ) , 

^-i = y; (22) 

where we have employed the identities H\(z) = z and Zerfc(z) = — zerfc(z) + 
(tt)- 1 ^- 2 . 

b. matching inner and outer solutions 
The inner solution 

p\x, y) = X 1 {x, y)Y 1 {x, y) + X^{x, y)Y^{x, y), (23) 
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in the limit of large \x\ and \y\, tends to the following functional form: 

For y > 0: 2{B 2 - D x )xy if x < 0, 2{B 2 + C x )xy if x > 0; 

For y < 0: 2(-A 2 - D x )xy if x < 0, 2(-A 2 + C x )xy if a; > 0. (24) 

We see that this outer limit of p 1 is proportional to xy in each of the four 
quadrants, The coefficient of proportionality may be specified independently in 
three of the four quadrants by choosing a suitable combination of A 2 , B 2 , C\, and 
D\. The fluid velocity at the outer scale may have two independently-specified 
discontinuities along the positive and negative x and y axes. The third indepen- 
dent discontinuity depends on the other two (the sum of all the discontinuities 
must of course equal zero). Since the integral condition (piTTj) may specify four 
independent proportionality conditions, and thus three independent velocity dis- 
continuities, we do not yet have a complete agreement between the inner and outer 
solutions. 

5 Conclusion 

The fact that one may have, in the outer limit, three independently-specified dis- 
continuities along the positive and negative x and y axes in the vicinity of the 
bathymetric saddle point, means that it is possible to specify different integral 
constraints for the circulation in three of the four quadrants. Extending this result 
to the larger-scale flow pattern in a basin, this means that flow velocity discon- 
tinuities predicted from (J§HTTj) across separatrices in the isobaths are not quite 
resolvable at an inner scale by solutions of the form given in (j22fJ2*Hj) in terms of 
the current analysis. Further investigation is required in order to obtain complete 
agreement between the solutions at the inner and outer scales. 

The solutions calculated in this paper are valid in the vicinity of saddle points 
which have separatrices crossing at right angles. For more general types of saddle 
point, the homogeneous version of (jHJ) is no longer separable: however, I anticipate 
that the matching of inner and outer solutions will eventually allow for the appro- 
priate specification of outer-scale flow discontinuities across separatrices. Hence 
the integral-constraint specification of the steady circulation proposed by Nilsson 
et al. j2] should be valid even for bathymetries containing saddle points. 
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